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1 

HYBRIDIZATION AND SEQUENCING OF NUCLEIC ACIDS 

GOVERNMENT RIGHTS 
The invention described herein arose in the course of or 
5 under Contract No. DE-FG03-92ER81275 (Grant No. 21012-92-11) between 

the Department of Energy and Affymax; and in the course of or under 
NIH Contract No. 1R01HG00813-01 . 

BACKGROUND OF THE INVENTION 

10 The present invention relates to the field of nucleic 

acid analysis, detection, and sequencing. More specifically, in 
one embodiment the invention provides improved techniques for 
synthesizing arrays of nucleic acids, hybridizing nucleic acids, 
detecting mismatches in a double- stranded nucleic acid composed of a 

15 single-stranded probe and a target nucleic acid, and determining the 

sequence of DNA or RNA or other polymers . 

It is important in many fields to determine the sequence 
of nucleic acids because, for example, nucleic acids encode the 
enzymes, structural proteins, and other effectors of biological 

20 functions. In addition to segments of nucleic acids that encode 

polypeptides, there are many nucleic acid sequences involved in 
control and regulation of gene expression. 

The human genome project is one example of a project 
using nucleic acid sequencing techniques. This project is directed 

25 toward determining the complete sequence of the genome of the human 

organism. Although such a sequence would not necessarily correspond 
to the sequence of any specific individual, it will provide 
significant information as to the general organization and specific 
sequences contained within genomic segments from particular 

30 individuals. The human genome project will also provide mapping 

information useful for further detailed studies. 

The need for highly rapid, accurate, and inexpensive 
sequencing technology is nowhere more apparent than in a demanding 
sequencing project. such as the human genome project. To complete the 

35 sequencing of a human genome will require the determination of 

approximately 3xl0 9 , or 3 billion, base pairs. 

The procedures typically used today for sequencing 
include the methods described in Sanger et al. ( Proc. Natl. Acad. 
Sci . USA (1977) 74:5463-5467, and Maxam et al., Methods in Enzymolocrv 

40 (1980) 65:499-559. The Sanger method utilizes enzymatic elongation 

with chain terminating dideoxy nucleotides. The Maxam and Gilbert 
method uses chemical reactions exhibiting base-specific cleavage 
reactions. Both methods require a large number of complex 
manipulations, such as isolation of homogeneous DNA fragments, 
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elaborate and tedious preparation of samples, preparation of a 
separating gel, application of samples to the gel, electrophoresing 
the samples on the gel, working up of the finished gel, and analysis 
of the results of the procedure. 
5 Alternative techniques have been proposed for sequencing 

a nucleic acid. PCT patent Publication No. 92/10588, incorporated 
herein by reference for all purposes, describes one improved 
technique in which the sequence of a labeled, target nucleic acid is 
determined by hybridization to an array of nucleic acid probes on a 

10 substrate. Each probe is located at a positionally distinguishable 

location on the substrate. When the labeled target is exposed to the 
substrate, it binds at locations that contain complementary 
nucleotide sequences. Through knowledge of the sequence of the 
probes at the binding locations, one can determine the nucleotide 

15 sequence of the target nucleic acid. The technique is particularly 

efficient when very large arrays of nuleic acid probes are utilized. 
Such arrays can be formed according to the techniques described in 
U.S. Patent No. 5,143,854 issued to Pirrung et al. See also U.S. 
application Serial No. 07/805,727, both incorporated herein by 

20 reference for all purposes. 

When the nucleic acid probes are of a length shorter than 
the target, one can employ a reconstruction technique to determine 
the sequence of the larger target based on affinity data from the 
shorter probes. See U.S. Patent No. 5,202,231 to Drmanac et al., 

25 and PCT patent Publication No. 89/10977 to Southern. One technique 

for overcoming this difficulty has been termed sequencing by 
hybridization or SBH. For example, assume that a 12-mer target DNA 
5 ' -AGCCTAGCTGAA is mixed with an array of all octanucleotide probes. 
If the target binds only to those probes having an exactly 

30 complementary nucleotide sequence, only five of the 65,536 octamer 

probes (3 ' -TCGGATCG, CGGATCGA, GGATCGAC, GATCGACT, and ATCGACTT) will 
hybridize to the target. Alignment of the overlapping sequences from 
the hybridizing probes reconstructs the complement of the original 
12-mer target: 

35 

TCGGATCG 
CGGATCGA 
GGATCGAC 
GATCGACT 
40 ATCGACTT 

TCGGATCGACTT 

While meeting with much optimism, prior techniques have 
also met with certain limitations. For example, practitioners have 
45 encountered substantial difficulty in analyzing probe arrays 
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hybridized to a target nucleic acid due to the hybridization of 
partially mismatched sequences, among other difficulties. The 
present invention provides significant advances in sequencing with 
such arrays . 

5 

SUMMARY OF THE INVENTION 
Improved techniques for synthesizing, hybridizing, 
analyzing, and sequencing nucleic acids (oligonucleotides) are 
provided by the present invention. 

10 According to one embodiment of the invention, a target 

oligonucleotide is exposed to a large number of immobilized probes 
of shorter length. The probes are collectively referred to as an 
"array." In the method, one identifies whether a target nucleic acid 
is complementary to a probe in the array by identifying first a core 

15 probe having high affinity to the target, and then evaluating the 

binding characteristics of all probes with a single base mismatch as 
compared to the core probe. If the single base mismatch probes 
exhibit a characteristic binding or affinity pattern, then the core 
probe is exactly complementary to at least a portion of the target 

20 nucleic acid. 

The method can be extended to sequence a target nucleic 
acid larger than any probe in the array by evaluating the binding 
affinity of probes that can be termed "left" and "right" extensions 
of the core probe. The correct left and right extensions of the 

25 core are those that exhibit the strongest binding affinity and/or 

a specific hybridization pattern of single base mismatch probes. 
The binding affinity characteristics of single base mismatch probes 
follow a characteristic pattern in which probe/target complexes with 
mismatches on the 3' or 5' termini are more stable than probe/target 

30 complexes with internal mismatches. The process is then repeated to 

determine additional left and right extensions of the core probe to 
provide the sequence of a nucleic acid target . 

In some embodiments, such as in diagnostics, a target is 
expected to have a particular sequence. To determine if the target 

35 has the expected sequence, an array of probes is synthesized that 

includes a complementary probe and all or some subset of all single 
base mismatch probes. Through analysis of the hybridization pattern 
of the target to such probes, it can be determined if the target has 
the expected sequence and, if not, the sequence of the target may 

4 0 optionally be determined. 

Kits for analysis of nucleic acid targets are also 
provided by virtue of the present invention. According to one 
embodiment, a kit includes an array of nucleic acid probes. The 
probes may include a perfect complement to a target nucleic acid. 

4 5 The probes also include probes that are single base substitutions of 
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the perfect complement probe. The kit may include one or more of the 
A, C, T, G, and/or U substitutions of the perfect complement. Such 
kits will have a variety of uses, including analysis of targets for 
a particular genetic sequence, such as in analysis for genetic 
5 diseases. 

A further understanding of the nature and advantages of 
the inventions herein may be realized by reference to the remaining 
portions of the specification and the attached drawings. 

10 BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 illustrates light-directed synthesis of 
oligonucleotides. A surface (2) bearing photoprotected hydroxyls 
(OX) is illuminated through a photolithographic mask (M x ) generating 
free hydroxyls (OH) in the photodeprotected regions. The hydroxyl 

15 groups are then coupled to a 5 ' -photoprotected deoxynucleoside 

phosphoramidite (e.g., T-X) . A new mask <M 2 ) is used to illuminate 
a new pattern on the surface, and a second photoprotected 
phosphoramidite (e.g., C-X) is then coupled. Rounds of illumination 
and coupling are repeated until the desired set of oligonucleotide 

20 probes is obtained. A target <R) is exposed to the oligonucleotides, 

optionally with a label (*) . The location(s) where the target binds 
to the array is used to determine the sequence of the target; 

Fig. 2 illustrates hybridization and thermal dissociation 
of oligonucleotides, showing a fluorescence scan of a target nucleic 

25 acid (5 ' -GCGTAGGC- fluorescein) hybridized to an array of probes. The 

substrate surface was scanned with a Zeiss Axioscop 20 microscope 
using 488 nm argon ion laser excitation. The fluorescence emission 
above 520 nm was detected using a cooled photomultiplier (Hamamatsu 
934-02) operated in photon counting mode. The signal intensity is 

30 indicated on the scale shown to the right of the image. The 

temperature is indicated to the right of each panel in °C; 

Fig. 3 illustrates the sequence specificity of 
hybridization. (A) is an index of the probe composition at each 
synthesis site. 3'-CGCATCCG surface immobilized probe (referred 

35 to herein as S-3 ' -CGCATCCG) was synthesized in stripes 1, 3, and 5, 

and the probe S-3 ' -CGCTTCCG was synthesized in stripes 2, 4, and 6. 
(B) is a fluorescence image showing hybridization of the substrate 
with a target nucleic acid (10 nM 5' -GCGTAGGC- fluorescein) . 
Hybridization was performed in 6X SSPE, 0.1% Triton X-100 at 15°C for 

40 15 min . (C) is a fluorescence image showing hybridization with a 

second nucleic acid (10 nM 5 ' -GCGAAGGC) added to the hybridization 
solution of (B) . (D) is a fluorescence image showing hybridization 
results after (1) high temperature dissociation of f luoresceinated 
targets from (C) ; and (2) incubation of the substrate with a target 

45 nucleic acid (10 nM 5' -GCGAAGGC) at 15°C for 15 min. (E) is a 
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fluorescence image showing hybridization with a second nucleic acid 
(10 nM 5 ' -GCGTAGGC) added to the hybridization solution of (D) ,- 

Fig. 4 illustrates combinatorial synthesis of 4 4 
tetranucleotides . In round 1, one-fourth of the synthesis area is 
5 activated by illumination through mask 1 for coupling of the first 

MeNPoc-nucleoside (T in this case). In cycle 2 of round 1, mask 2 
activates a different one quarter section of the synthesis substrate, 
and a different nucleoside (C) is coupled. Further lithographic 
subdivisions of the array and chemical couplings generate the 
10 complete set of 256 tetranucleotides; 

Figs. 5A and 5B illustrate hybridization to an array 
of 256 octanucleotides . Fig. 5A is a fluorescence image following 
hybridization of the array with a target nucleic acid (10 nM 
5' -GCGGCGGC- fluorescein) in 6X SSPE, 0.1% Triton X-100 for 15 min. at 
15 15 °C. Fig. 5B is a matrix de -coder showing where each probe made 

during the synthesis of S-3 ' -CG (A+G+C+T) 4 CG is located. The site 
containing the probe sequence S-3 ' -CGCGCCCG is shown as a dark area. 
The combinatorial synthesis notation used herein is fully described 
in U.S. application Serial No. 07/624,120, incorporated herein by 
20 reference for all purposes.; 

Figs. 6A to 6C illustrate a technique for sequencing 
a n-mer target using k-mer probes. Fig. 6A illustrates a target 
hybridized to a probe on a substrate. Figs. 6B and 6C illustrate 
plots of normalized binding affinity vs. mismatch positiori; 
25 Fig. 7 illustrates a fluorescence image -of a 

hybridization experiment; 

Fig. 8 illustrates hybridization events graphically as a 
function of singe base mismatch; 

Fig. 9 illustrates fluorescence intensity as a function 
30 of pairs of mismatches; 

Fig. 10 illustrates a fluorescence image of a single base 
mismatch experiment; 

Figs. 11A to 11C illustrate various single base mismatch 

profiles ; 

35 Figs. 12A to 12D illustrate a process for determining the 

nucleotide sequence of an n-member (the number of monomers in the 
nucleotide) target oligonucleotide based on hybridization results 
from shorter k-member probes. In particular, Figs. 12A to 12D 
illustrate application of the present method to sequencing a 10-base 

4 0 target with 4 -base probes; 

Fig. 13 illustrates a computer system for determining 
nucleotide sequence; 

Fig. 14 illustrates a computer program for mismatch 
analysis and for determining the nucleotide sequence of a target 

45 nucleic acid; 
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Figs. 15A and 15B illustrate a computer program for 
determining the nucleotide sequence of a target nucleic acid by- 
selecting among several possibilities, and an example of a scoring 
routine for use in this computer program. 
5 Fig. 16 illustrates a directed graph for use in 

determining nucleotide sequence. 

Figs. 17A and 17B illustrate wild-type and mutation 
analysis using single base mismatch profiles ; 

Fig. 18 is a fluorescence image of a single base mismatch 

10 test; and 

Figs. 19A to 19D illustrate a technique for nucleic acid 
sequence identification. 



15 



DESCRIPTION OF THE PREFERRED EMBODIMENTS 
CONTENTS 



A. Synthesis 

B. Hybridization 

20 C. Mismatch Analysis 

D. Applications 

E. Conclusion 



25 Definitions 

Probe - A molecule of known composition or monomer 
sequence, typically formed on a solid surface, which is or may be 
exposed to a target molecule and examined to determine if the probe 
has hybridized to the target. A "core" probe is a probe that 

30 exhibits strong affinity for a target. An "extension" probe is a 

probe that includes all or a portion of a core probe sequence plus 
one or more possible extensions of the core probe sequence. The 
present application refers to "left 11 extensions as an extension at 
the 3' -end of a probe and a "right" extension refers to an extension 

35 at the 5' -end of a probe, although the opposite notation could 

obviously be adopted. 

Target - A molecule, typically of unknown composition or 
monomer sequence, for which it is desired to study the composition or 
monomer sequence. A target may be a part of a larger molecule, such 

40 as a few bases in a longer nucleic acid. 

n-Base Mismatch - A probe having n monomers therein that 
differ from the corresponding monomers in a core probe, wherein n is 
one or greater. 

A, T, C, G, U - Are abbreviations for the nucleotides 

45 adenine, thymine, cytosine, guanine, and uridine, respectively. 
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Library - A collection of nucleic acid probes of 
predefined nucleotide sequence, often formed in one or more 
substrates, which are used in hybridization studies of target nucleic 
acids . 

5 

A. Synthesis 

A method for a light-directed oligonucleotide synthesis 
is depicted in Fig. 1. Such strategies are described in greater 
detail in U.S. Patent No. 5,143,854, assigned to the assignee of the 
10 present inventions and incorporated herein by reference for all 

purposes . 

In the light-directed synthesis method illustrated in 
Fig. l, a surface 2 derivatized with a photolabile protecting group 
or groups (X) is illuminated through a photolithographic mask M 1 

15 exposing reactive hydroxyl {OH) groups. The first (T-X) of a 

series of phosphoramidite activated nucleosides (protected at the 
5' -hydroxyl with a photolabile protecting group) is then exposed to 
the entire surface. Coupling only occurs at the sites that were 
exposed to light during the preceding illumination. 

20 After the coupling reaction is complete, the substrate 

is rinsed, and the surface is again illuminated through a new or 
translated mask: M 2 to expose different groups for coupling. A new 
phosphoramidite activated nucleoside C-X (again protected at the 
5' -hydroxyl with a photolabile protecting group) is added and coupled 

25 to the exposed sites. The process is repeated through cycles of 

photodeprotection and coupling to produce a desired set of 
oligonucleotide probes on the substrate. Because photolithography 
is used, the process can be miniaturized. Furthermore, because 
reactions only occur at sites spatially addressed by light, the 

30 nucleotide sequence of the probe at each site is precisely known, and 

the interaction of oligonucleotide probes at each site with target 
molecules (either target nucleic acids or, in other embodiments, 
proteins such as receptors) can be assessed. 

Photoprotected deoxynucleosides have been developed for 

35 this process including 5 ' -0- (a-methyl-6-nitropiperonyloxycarbonyl) -N- 

acyl -2 ' -deoxynucleosides , or MeNPoc-N-acyl -deoxynucleosides, MeNPoc- 
dT, MeNPoc-dC it>u , MeNPoc-dG PAC , and MeNPoc-dA PAC . Protecting group 
chemistry is disclosed in greater detail in PCT patent Publication 
No. 92/10092 and U.S. application Serial Nos . 07/624,120, filed 

40 December 6, 1990, and 07/971,181, filed November 2, 1992, both 

assigned to the assignee of this invention and incorporated herein by 
reference for all purposes. 
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Examples 

1 . Protecting Groups 

Because the bases have strong 7r- tt* transitions in 
the 280 nm region, the deprotection wavelength of photoremovable 
5 protecting groups should be at wavelengths longer than 280 nm to 

avoid undesirable nucleoside photochemistry. In addition, the 
photodeprotection rates of the four deoxynucleosides should be 
similar, so that light will equally deprotect hydroxy Is (or other 
functional groups, such as sulfhydryl or amino groups) in all 
10 illuminated synthesis sites. 

To meet these criteria, a set of 5' -0- (af-methyl-6- 
nitropiperonyloxycarbonyl) -N-acyl-2' -deoxynucleosides (MeNPoc-N-acyl- 
deoxynucleosides) has been developed for light-directed synthesis, 
and the photokinetic behavior of the protected nucleosides has been 
15 measured. The synthetic pathway for preparing 5 ' -0' - (a-methyl-6- 

nitropiperonyloxycarbonyl) -N-acyl-2' -deoxynucleoside phosphoramidites 
is illustrated in Scheme I. 



O 




35 

Scheme I 



In the first step, an N-acyl-2' -deoxynucleoside was reacted with 
1- (2-nitro-4 , 5-methylenedioxyphenyl) -ethan- 1-chlorof ormate to yield 

40 5' -MeNPoc-N-acyl-2' -deoxynucleoside . In the second step, the 

3'-hydroxyl was reacted with 2-cyanoethyl-N,N' -diisopropylchloro- 
phosphoramidite using standard procedures to yield the 5'-MeNPoc-N- 
acyl-2 ' -deoxynucleoside -3 ' -O-diisopropylchlorophosphoramidites . 
These reagents were stable for long periods when stored dry under 

45 argon at 4°C. 
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A 0.1 mM solution of each of the four deoxynucleosides, 
MeNPoc-dT, MeNPoc-dC ibu , MeNPoc-dG PAC , and MeNPoc-dA PAC was prepared 
in dioxane. Aliquots (200 mL) were irradiated with 14.5 mW/cm 2 365 
nm light in a narrow path (2 mm) quartz cuvette for various times. 
5 Four or five time points were collected for each base, and the 

solutions were analyzed for loss of starting material with an HPLC 
system at 280 nm and a nucleosil 5-Cg HPLC column, eluting with a 
mobile phase of 60% (v/v) in water containing 0.1% (v/v) TFA 
(MeNPoc-dT required a mobile phase of 70% (v/v) methanol in water) . 

10 Peak areas of the residual MeNPoc -N-acyl-deoxynucleoside were 

calculated, yielding photolysis half-times of 28 s, 31 s, 27 s, and 
18 s for MeNPoc-dT, MeNPoc-dC ibu , MeNPoc-dG PAC , and MeNPoc-dA PAC , 
respectively. In subsequent lithographic experiments, illumination 
times of 4.5 min. (9*t 1 / 2 MeNPoc_dC ) led to more than 99% removal of 

15 MeNPoc protecting groups. 

In a light-directed synthesis, the overall synthesis 
yield depends on the photodeprotection yield, the photodeprotection 
contrast, and the chemical coupling efficiency. Photokinetic 
conditions are preferably chosen to ensure that photodeprotection 

20 yields are over 99%. Unwanted photolysis in normally dark regions of 

the substrate can adversely affect the synthesis fidelity but can be 
minimized by using lithographic masks with a high optical density 
(5 ODU) and by careful index matching of the optical surfaces. 
Condensation efficiencies of DMT-N-acyl-deoxynucleoside 

25 phosphoramidites to the glass substrates have been measured in the 

range of 95% to 99%. The condensation efficiencies of the MeNPoc-N- 
acyl-deoxynucleoside phosphoramidites have also been measured at 
greater than 9 0%, although the efficiencies can vary from synthesis 
to synthesis and should be monitored. 

30 

2 . Coupling Efficiency Measurements 

To investigate the coupling efficiencies 
of the photoprotected nucleosides, each of the four MeNPoc-amidites 
was first coupled to a substrate (via DMT chemistry) . A region of 

35 the substrate was illuminated, and a MeNPoc -phosphoramidite was 

added without a protective group. A new region of the substrate was 
then illuminated; a fluorescent deoxynucleoside phosphoramidite 
(FAM-phosphoramidite Applied Biosystems) was coupled; and the 
substrate was scanned for signal. If the f luorescently labeled 

4 0 phosphoramidite reacts at both the newly exposed hydroxyl groups 

and the previously unreacted hydroxyl groups, then the ratio of 
fluorescence intensities between the two sites provides a measure 
of the coupling efficiency. This measurement assumes that surface 
photolysis yields are near unity. The chemical coupling yields using 

45 this or similar assays are variable but high, ranging between 80-95%. 
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10 



10 

in a separate assay, the chemical coupling efficiencies 
were measured on hexaethyleneglycol derivatized substrate. First, 
a glycol linker was detritylated and a MeNPoc-deoxynucleoside-O- 
cyanoethylphosphoramidite coupled to the resin without capping. 
Next, a DMT-deoxyucleoside-cyanoethylphosphoramidite (reporter- 
amidite) was coupled to the resin. The reporter-amidite couples 
to any unreacted hydroxyl groups from the first step. The trityl 
effluents were collected and quantified by absorption spectroscopy. 
Effluents were also collected from the lines immediately after the 
MeNPoc-phosphoramidite coupling to measure residual trityl left in 
the delivery lines. In this assay, the coupling efficiencies are 
measured assuming a 100% coupling efficiency of the reporter-amidite. 
The coupling efficiencies of the MeNPoc-deoxyribonucleoside-O- 
cyanoethylphosphoramidites to the hexaethyleneglycol linker and the 
15 efficiencies of the sixteen dinucleotides were measured and were 

indistinguishable from DMT- deoxynucleo side phosphoramidites . 

3 - Spatiall y Directed Synthesis of an Oligonucleotide Probe 
To initiate the synthesis of an oligonucleotide probe, 

20 substrates were prepared, and MeNPoc-dC ibu -3 ' -O-phosphoramidite was 

attached to a synthesis support through a synthetic linker. Regions 
of the support were activated for synthesis by illumination through 
800 x 1280 jim apertures of a photolithographic mask. Seven 
additional phosphoramidite synthesis cycles were performed (with 

25 the corresponding DMT protected deoxynucleosides) to generate the 

S-3' -CGCATCCG. Following removal of the phosphate and exocyclic 
amine protecting groups with concentrated NH 4 0H for 4 hours at 
room temperature, the substrate was mounted in a water jacketed 
thermostatically controlled hybridization chamber. This substrate 

30 was used in the mismatch experiments referred to below. 

B . Hybridization 

Oligonucleotide arrays can be used in a wide variety of 
applications, including hybridization studies. In a hybridization 

35 study, the array can be exposed to a receptor (R) of interest, as 

shown in Fig 1. The receptor can be labelled with an appropriate 
label (*), such as fluorescein. The locations on the substrate where 
the receptor has bound are determined and, through knowledge of the 
sequence of the oligonucleotide probe at that location one can then 

40 determine, if the receptor is an oligonucleotide, the sequence of the 

receptor. 

Sequencing by hybridization (SBH) is most efficiently 
practiced by attaching many probes to a surface to form an array in 
which the identity of the probe at each site is known. A labeled 
45 target DNA or RNA is then hybridized to the array, and the 
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11 

hybridization pattern is examined to determine the identity of all 
complementary probes in the array. Contrary to the teachings of the 
prior art, which teaches that mismatched probe/target complexes are 
not of interest, the present invention provides an analytical method 
5 in which the hybridization signal of mismatched probe/target 

complexes identifies or confirms the identity of the perfectly 
matched probe/target complexes on the array. 

Arrays of oligonucleotides are efficiently generated for 
the hybridization studies using light-directed synthesis techniques. 

10 As discussed below, an array of all tetranucleotides was produced in 

sixteen cycles, which required only 4 hours to complete. Because 
combinatorial strategies are used, the number of different compounds 
on the array increases exponentially during synthesis, while the 
number of chemical coupling cycles increases only linearly. For 

15 example, expanding the synthesis to the complete set of 4^ (65,536) 

octanucleotides adds only 4 hours (or less) to the synthesis due to 
the 16 additional cycles required. Furthermore, combinatorial 
synthesis strategies can be implemented to generate arrays of any 
desired probe composition. For example, because the entire set of 

20 dodecamers (4 12 ) can be produced in 48 photolysis and coupling cycles 

or less (b n compounds requires no more than b x n cycles) , any subset 
of the dodecamers (including any subset of shorter oligonucleotides) 
can be constructed in 48 or fewer chemical coupling steps. The 
number of compounds in an array is limited only by the density of 

25 synthesis sites and the overall array size. The present invention" 

has been practiced with arrays with probes synthesized in square 
sites 25 microns on a side. At this resolution, the entire set of 
65,536 octanucleotides can be placed in an array measuring only 
0.64 cm 2 . The set of 1,048,576 dodecanucleotides requires only a 

30 2.56 cm 2 array at this individual probe site size. 

The success of genome sequencing projects depends on 
efficient DNA sequencing technologies. Current methods are highly 
reliant on complex procedures and require substantial manual effort. 
SBH offers the potential for automating many of the manual efforts in 

35 current practice. Light-directed sytnhesis offers an efficient means 

for large scale production of miniaturized arrays not only for SBH 
but for many other applications as well. 

Although oligonucleotide arrays can be used for primary 
sequencing applications, many diagnostic methods involve the analysis 

40 of only a few nucleotide positions in a target nucleic acid sequence. 

Because single base changes cause multiple changes in the 
hybridization pattern of the target on a probe array, the 
oligonucleotide arrays and methods of the present invention enable 
one to check the accuracy of previously elucidated DNA sequences, or 

45 to scan for changes or mutations in certain specific sequences within 
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a target nucleic acid. The latter as is important, for example, for 
genetic, disease, quality control, and forensic analysis. With an 
octanucleotide probe set, a single base change, in a target nucleic 
acid can be detected by the loss of eight perfect hybrids, and the 
5 generation of eight new perfect hybrids. The single base change can 

also be detected through altered mismatch probe/target complex 
formation on the array. Perhaps even more surprisingly, such single 
base changes in a complex nucleic acid dramatically alter the overall 
hybridization pattern of the target to the array. According to the 

10 present invention such changes in the overall hybridization pattern 

are used to actually simplify the analysis. 

The high information content of light-directed 
oligonucleotide arrays greatly benefits genetic diagnostic testing. 
Sequence comparisons of hundreds to thousands of different mutations 

15 can be assayed simultaneously instead of in a one-at-a-time format. 

Arrays can also be constructed to contain genetic markers for the 
rapid identification of a wide variety of pathogenic organisms, 
and to study the sequence specificity of RNA/RNA, RNA/DNA, 
protein/RNA or protein/DNA, interactions. One can use non Watson- 

20 Crick oligonucleotides and novel synthetic nucleoside analogs for 

antisense, triple helix, or other applications. Suitably protected 
RNA monomers can be employed for RNA synthesis, and a wide variety of 
synthetic and non-naturally occurring nucleic acid analogues can be 
used, depending upon the motivations of the practitioner. See, e.g., 

25 PCT patent Publication Nos . 91/19813, 92/05285, and 92/14843, 

incorporated herein by reference. In addition, the oligonucleotide 
arrays can be used to deduce thermodynamic and kinetic rules 
governing the formation and stability of oligonucleotide complexes. 

30 Examples 

1 . Hybridization of Targets to Surface Oligonucleotides 
The support bound octanucleotide probes discussed 
above were hybridized to a target of 5 'GCGTAGGC- fluorescein in 
the hybridization chamber by incubation for 15 minutes at 15°C. 

35 The array surface was then interrogated with an epif luorescence 

microscope (488 nm argon ion excitation) . The fluorescence image of 
this scan is shown in Fig. 2. The fluorescence intensity pattern 
matches the 800 X 1280 /im stripe used to direct the synthesis of the 
probe. Furthermore, the signal intensities are high (four times over 

40 the background of the glass substrate) , demonstrating specific 

binding of the target to the probe. 

The behavior of the target-probe complex was investigated 
by increasing the temperature of the hybridization solution. After a 
10 minute equilibration at each temperature, the substrate was 

45 scanned for signal. The duplex melted in the temperature range 
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expected for the sequence under study (T m ~28°C obtained from the rule 
T m = [2° (A+T)+4° (G+C) ] ) . The probes in the array were stable to 
temperature denaturation of the target-probe complex as demonstrated 
by rehybridization of target DNA. 

5 

2 . Sequence Specificity of Target Hybridization 

To demonstrate the sequence specificity of target 
hybridization, two different probes were synthesized in 800 x 1280 fim 
stripes. Fig. 3A identifies the location of the two probes. The 

0 probe S - 3 ' - CGCATCCG was synthesized in stripes l, 3 and 5. The probe 

S-3' -CGCTTCCG was synthesized in stripes 2, 4 and 6. Fig. 3B shows 
the results of hybridizing a 5 ' -GCGTAGGC- fluorescein target to the 
substrate at 15°C. Although the probes differ by only one internal 
base, the target hybridizes specifically to its complementary 

5 sequence (-500 counts above background in stripes 1, 3 and 5) with 

little or no detectable signal in positions 2, 4 and 6 (-10 counts). 
Fig. 3C shows the results of hybridization with targets to both 
sequences. The signal in all positions in Fig. 3C illustrates that 
the absence of signal in Fig. 3B is due solely to the instability 

0 of the single base mismatch. Although the targets are present in 

equimolar concentrations, the ratio of signals in stripes 2, 4 and 6 
in Fig. 3B are approximately 1.6 times higher than the signals in 
regions l, 3 and 5. This duplex has a slightly higher predicted T m 
than the duplex comprising regions 2, 4 and 6. The duplexes were 

5 dissociated by raising the temperature to 45°C for 15 minutes, and 

the hybridizations were repeated in the reverse order (Figs. 3D and 
3E) , demonstrating specificity of hybridization in the reverse 
direction . 



3 • Combinatorial Synthesis of, and Hybridization 
of a Nucle ic Acid Target to, a Probe Matrix 

In a light-directed synthesis, the location 

and composition of products depends on the pattern of illumination 

and the order of chemical coupling reagents (see Fodor et al., 

Science (1991) 251:767-773, for a complete description). Consider 

the synthesis of 256 tetranucleotides , as illustrated in Fig. 4. 

Mask 1 activates one fourth of the substrate surface for coupling 

with the first of four nucleosides in the first round of synthesis. 

In cycle 2, mask 2 activates a different quarter of the substrate for 

coupling with the second nucleoside. The process is continued to 

build four regions of mononucleotides. The masks of round 2 are 

perpendicular to those of round 1, and each cycle of round 2 

generates four new dinucleotides . The process continues through 

round 2 to form sixteen dinucleotides as illustrated in Fig. 4. The 

masks of round 3 further subdivide the synthesis regions so that each 
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coupling cycle generates 16 trimers . The subdivision of the 
substrate is continued through round 4 to form the tetranucleot ides . 
The synthesis of this probe matrix can be compactly represented in 
polynomial notation as (A+OG+T) 4 . Expansion of this polynomial 
5 yields the 256 tetranucleotides . 

The application of an array of 256 probes synthesized by 
light-directed combinatorial synthesis to generate a probe matrix is 
illustrated in Fig. 5A. The polynomial for this synthesis is given 
by: 3' -CG(A+G+C+T) 4 CG. The synthesis map is given in Fig. 5B . All 

10 possible tetranucleotides were synthesized flanked by CG at the 3'- 

and 5 '-ends. Hybridization of target 5' -GCGGCGGC- fluorescein to this 
array at 15°C correctly yielded the S-3 ' -CGCCGCCG complementary probe 
as the most intense position (2,698 counts). Significant intensity 
was also observed for the following mismatches: S-3 ' -CGCAGCCG 

15 (554 counts), S - 3 ' - CGCCGACG (317 counts), S - 3 ' - CGCCGTCG (272 counts), 

S-3' -CGACGCCG (242 counts), S-3 ' -CGTCGCCG (203 counts), S - 3 ' - CGCCCCCG 
(180 counts), S-3' -CGCTGCCG (163 counts), S - 3 ' - CGCCACCG (125 counts), 
and S-3' -CGCCTCCG (78 counts). 



20 C. Mismatch Analysis 

The arrays discussed above can be utilized in the present 
method to determine the nucleic acid sequence of an oligonucleotide 
of length n using an array of probes of shorter length k. Fig. 6 
illustrates a simple example. The target has a sequence 5'-XXYXY-3', 

25 where X and Y are complementary nucleic acids such as A and T or C 

and G. For discussion purposes, the illustration in Fig. 6 is 
simplified by using only two bases and very short sequences, but the 
technique can easily be extended to larger nucleic acids with, for 
example, all 4 RNA or DNA bases. 

30 The sequence of the target is, generally, not known 

ab initio . One can determine the sequence of the target using the 
present method with an array of shorter probes. In this example, an 
array of all possible X and Y 4-mers is synthesized and then used to 
determine the sequence of a 5-mer target. 

35 Initially, a "core" probe is identified. The core probe 

is exactly complementary to a sequence in the target using the 
mismatch analysis method of the present invention. The core probe is 
identified using one or both of the following criteria: 

1 . The core probe exhibits stronger binding 

40 affinity to the target than other probes, typically the 

strongest binding affinity of any probe in the array 
(that has not been identified as a core probe in a 
previous cycle of analysis) . 

2. Probes that are mismatched with the target, 
45 as compared to the core probe sequence, exhibit a 
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characteristic pattern, discussed in greater detail 
below, in which probes that mismatch at the 3'- and 
^ 5' -end of the probe bind more strongly to the target than 
probes that mismatch at interior positions. 
5 In this particular example, selection criteria #1 identifies a core 

4-mer probe with the strongest binding affinity to the target that 
has the sequence 3'-YYXY, as shown in Fig. 6A, where the probe is 
illustrated as having hybridized to the target. The probe 3'-YYXY 
(corresponding to the 5'-XXYX position of the target) is, therefore, 

10 chosen as the "core" probe. 

Selection criteria #2 is utilized as a "check" to ensure 
the core probe is exactly complementary to the target nucleic acid. 
The second selection criteria evaluates hybridization data (such as 
the fluorescence intensity of a labeled target hybridized to an array 

15 of probes on a substrate, although other techniques are well known to 

those of skill in the art) of probes that have single base mismatches 
as compared to the core probe. In this particular case, the core 
probe has been selected as S-3'-YYXY. The single base mismatched 
probes of this core probe are: S-3'-XYXY, S-3'-YXXY, S-3'-YYYY, and 

20 S-3'-YYXX. The binding affinity characteristics of these single base 

mismatches are utilized to ensure that a "correct" core has been 
selected, or to select the core probe from among a set of probes 
exhibiting similar binding affinities. 

An illustrative, hypothetical plot of expected binding 

25 affintity versus mismatch position is provided in Fig. 6B. The 

binding affinity values (typically fluorescence intensity of labeled 
target hybridized to probe, although many other factors relating to 
affinity may be utilized) are all normalized to the binding affinity 
of S-3'-YYXY to the target, which is plotted as a value of l on the 

30 left hand portion of the graph. Because only two nucleotides are 

involved in this example, the value plotted for a probe mismatched at 
position 1 (the nucleotide at the 3' -end of the probe) is the 
normalized binding affinity of S-3'-XYXY. The value plotted for 
mismatch at position 2 is the normalized affinity of S-3'-YXXY. The 

35 value plotted for mismatch at position 3 is the normalized affinity 

of S-3'-YYYY, and the value plotted for mismatch position 4 is the 
normalized affinity of S-3'-YYXX. As noted above, "affinity" may be 
measured in a number of ways including, for example, the number of 
photon counts from fluorescence markers on the target. 

4 0 The affinity of all three mismatches is lower than the 

core in this illustration. Moreover, the affinity plot shows that a 
mismatch at the 3 '-end of the probe has less impact than a mismatch 
at the 5' -end of the probe in this particular case, although this may 
not always be the case. Further, mismatches at the end of the probe 

45 result in less disturbance than mismatches at the center of the 
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10 



probe. These features, which result in a "smile" shaped graph when 
plotted as shown in Fig. 6B, will be found in most plots of single 
base mismatch after selection of a "correct" core probe, or after 
accounting for a mismatched probe that is a core probe with respect 
to another portion of the target sequence. This information will be 
utilized in either selecting the core probe initially or in checking 
to ensure that an exactly matched core probe has been selected. 
Of course, in certain situations, as noted in Section B above, 
identification of a core is all that is required such as in, 
for example, forensic or genetic studies, and the like. 

In sequencing studies, this process is then repeated for 
left and/or right extensions of the core probe. In the example 
illustrated in Fig. 6, only right extensions of the core probe are 
possible. The possible 4-mer extension probes of the core probe 
15 are 3'-YXYY and 3'-YXYX. Again, the same selection criteria are 

utilized. Between 3'-YXYY and 3'-YXYX, it would normally be found 
that 3'-YXYX would have the strongest binding affinity, and this 
probe is selected as the correct probe extension. This selection may 
be confirmed by again plotting the normalized binding affinity of 
20 probes with single base mismatches as compared to the core probe. 

A hypothetical plot is illustrated in Fig. 6C. Again, the 
characteristic "smile" pattern is observed, indicating that the 
"correct" extension has been selected, i.e., 3'-YXYX. From this 
information, one would correctly conclude that the sequence of the 
25 target is 5'-XXYXY. 



30 



Examples 

1 - (A+T) e A rray and Single Base Mismatch Stabilities 
A 20-step, 4-replica combinatorial synthesis was 
performed using MenPoc-dA and MenPoc-dT. The lithographic masks were 
chosen such that each member of a set of 256 octanucleot ides was 
synthesized in four separate locations on the 1.28 x 1.28 cm array, 
yielding 1024 different synthesis sites, each containing an 
octanucleotide probe, each site 400 x 400 in size. Following 
35 synthesis and phenoxyacetyl deprotection of the dA amine, the 

substrate was mounted in a thermostatically regulated staining and 
flow cell, incubated with 1 nM 5 ' -AAAAAAAA- fluoresce in at 15°C, and 
then scanned in a Zeiss epif luorescence microscope. The resulting 
fluorescent image is shown in Fig. 7. 
40 Fluorescence intensities of the hybridization events as a 

function of single base mismatch are provided graphically in Fig. 8. 
Each of the four independent intensities for each octanucleotide 
probe that differs from the core probe at a single base is plotted. 
Position zero mismatch {i.e., the perfect complement 3 ' -TTTTTTTT) i 
45 the brightest position on the array at -9 00 counts; the background 



is 
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signal of this array is approximately 220 counts. Mismatch position 
1 (at the 3' -end of the probe) is the next brightest at -760 counts. 
A "smile" or "U" shaped curve of the following, positions indicates 
the relative stability of the mismatches at each position of the 
5 probe/target complex. This "mismatch family" characteristizes 

nucleic acid interaction with an array of probes and provides or 
confirms the identification of the target sequence. The mismatches 
at positions 3, 4, 5 and 6 are more destabilizing and yield 
intensities virtually indistinguishable from background. 

10 The mismatch at position l (the point where the 3 '-end of the 

octanucleotide is tethered to the substrate) is less destablizing 
than the corresponding mismatch at position 8 (the free 5' -end). 
The uniformity of the array synthesis and the target hybridization 
is reflected in the low variation of intensities between the four 

15 duplicate synthesis sites. 

The method of the present invention can also utilize 
information from target hybridization to probes with two or more 
mismatches. Fluorescence intensities as a function of pairs of 
mismatches are presented in Fig. 9. In this case, the intensity data 

20 have been normalized so that a perfect match has intensity l. For 

example, the data at index 1,8 corresponds to mismatches at each end 
of the probe/target duplex. The diagonal (index 1,1 to 8,8) 
corresponds to the single mismatches illustrated in Fig. 8. The 
highest intensities correspond to single and pairs of mismatches at 

25 the ends of the probe/target complex. 

2 . (G+T) 8 Array and Sequence Reconstruction 

An octanucleotide array of MenPoc-dG and MenPoc-dT 
was synthesiszed. The format of the synthesis was similar to 

30 that for the (A+T) 8 array, discussed above, and resulted in 

256 octanucleotides of G and T in replicates of four (1024 total) . 
After final deprotection and attachment to a temperature-controlled 
(15°C) hybridization chamber, the probe array was incubated with 
1 nM 5' -AACCCAAACCC-f luorescein target and scanned. The resulting 

35 image is given in Fig. 10. Four distinct but overlapping, perfectly 

complementary octanucleotide hybridizations are expected: 
3 ' -TTGGGTTT, TGGGTTTG, GGGTTTGG, and GGTTTGGG . As shown herein, the 
moderate stablility of probe/target complexes with single base pair 
mismatches generates families of probes with moderate signals. 

40 A cursory inspection of the many intense features of Fig. 10 revealed 

a complex pattern. 

The reconstruction heuristic provided by the present 
invention effectively utilizes the complex data pattern in Fig. 10. 
The algorithm assumes as a general rule that perfectly matched 

45 probe/target complexes have higher fluorescence intensities, and 
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perfect matches and related single base mismatch typically form a 
profile similar to that shown in Fig. 6. 

The probe with the highest intensity should be a perfect 
match to the target. Corresponding mismatch profiles are shown in 
5 Figs. 11A to 11C. One first plots the mismatch profile for the probe 

with the highest intensity <S-3 ' -TGGGTTTG in this case) to verify 
that the probe is exactly complementary to the target. Assuming that 
this probe is complementary to a fragment of the target, we consider 
"extending" a base on the 3 '-end of the target. In this case, there 

10 are two probe choices. One of the two 8-mer probes S - 3 ' - GGGTTTGT and 

S-3 ' -GGGTTTGG, will be exactly complementary to the target nucleic 
acid. The mismatch profile for each of these two probes, as well as 
for probe S-3 ' -TGGGTTTG, is shown with intensity values in Fig. 11A. 
Note that the probe S-3 ' -GGGTTTGG has the mismatch profile most 

15 similar to that of probe S-3 ' -TGGGTTTG (a typical "smile" plot). 

Therefore, one will conclude that the correct extension probe is 
S-3' -GGGTTTGG. 

Fig. 11B shows repetition of this process to evaluate the 
3 '-end of the target sequence. Because the probe S- 3 ' -GGTTTGG has a 

20 smile-shaped mismatch profile most like the core S-3 ' -GGGTTTGG, and 

because the probe S-3 ' -GGTTTGGT does not, one will correctly conclude 
that the probe S- 3 ' -GGGTTTGG is the correct extension probe. This 
process can be repeated until neither profile has the correct shape, 
or the absolute intensity is well below that of the highest 

25 intensity, indicating that the "end" of the target has been reached. 

A similar method provides the sequence of the target extending to the 
5' -end. Fig. 11C shows the mismatch curves for all the perfectly 
matched probes; each curve has the consistent shape predicted for 
this target . 

30 The techniques described above can of course be readily 

extended to nucleic acids of any length, as illustrated in the 
various panels of Figs. 12A to 12D. As shown in Fig. 12A, a 10-mer 
target is to be sequenced, and the sequence is indicated by 
5 ' ~ N i N 2 N 3 N 4 N 5 N 6 N 7 N 8 N 9 N l0" 3 ' ' wnere N is" any nucleotide or nucleic 

35 acid monomer, and the subscript indicates the nucleotide position in 

the probe, with 3 indicating the 3 '-end terminal monomer. Those of 
skill recognize that, if the probes were synthesized with the 5' -end 
attached to the substrate, the method of the invention can be applied 
with appropriate modification. 

4 0 An array of shorter oligonucleotides can be used to 

sequence a larger nucleotide according to one aspect of the present 
invention. In the particular example shown in Figs. 12A to 12D, 4- 
mers (oligonucleotide probes 4 monomers in length) are used to 
sequence the unknown 10-mer target. In practice, longer probes and 

45 targets will typically be employed, but this illustrative example 
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facilitates understanding of the invention. A single member of the 

4- mer array is shown in Fig. 12A and has the sequence S-3 ' -P 3 P 4 P 5 P 6 , 
where the various P (probe) nucleotides will be selected from the 
group of A, T, C ( U, G, and other monomers, depending on the 

5 application, and the subscript indicates position relative to the 

target. For discussion purposes, the hybridization data are presumed 
to be available from a single array. However, one can utilize 
multiple arrays, arrays synthesized at different times, or even 
individual probes to practice the method. As noted, the probe length 

10 of 4 is selected to facilitate discussion; in practice, longer probes 

will typically be employed. 

S-3' -P 3 P 4 P 5 P 6 is selected as a core probe from the array 
due to its exhibition of a strong binding affinity to the target and 
a correct mismatch profile. In the array of all 4-mers, the sequence 

15 S-3 ' - p 3 p 4 p 5p6 is chosen as the core sequence, because when a 

f luorescein-labeled target (shown as 5' -N 1 N 2 N 3 N 4 N 5 N 6 N 7 N 8 N 9 N l0 *-3' in 
Fig. 12A) is exposed to the substrate, the target hybridizes to the 
probe, as indicated by the arrows in Fig. 12A, and high fluorescence 
intensity (i.e., a large number of photon counts) is observed in the 

20 portion of the substrate containing the probe S-3 ' -P 3 P 4 P 5 P 6 , as 

compared to other portions of the substrate. Normally, the sequence 
exhibiting the strongest binding affinity will be chosen as the first 
core sequence . 

One preferably verifies whether the first selected core 

25 sequence is a perfect- complement to the target by examining the 

fluorescence intensity of probes in the array that differ from the 
core probe at a single base. Fig. 12B qualitatively illustrates a 
typical plot of relative intensity of single base mismatches versus 
position of the mismatch for the S-3 ' -P 3 P 4 P 5 P 6 core probe. As a 

30 simple example, assume that, in the sequence S-3 ' -P 3 P 4 P 5 P 6 , the 

nucleotide C is not present. Fig. 12B illustrates in a qualitative 
way the normalized fluorescence intensity of probes that differ from 
the core sequence probe by substitution of C into the sequence 

5- 3' -P 3 P 4 P 5 P 6 and in which none of the C-containing mismatched probes 
35 is exactly complementary to another sequence in the target. 

Accordingly, Fig. 12B plots the relative fluorescence intensity of 
the probe set : 

S-3'-CP 4 P 5 P 6 , 
S-3'-P 3 CP 5 P 6 , 
40 S-3' -P 3 P 4 CP 6 , and 

S-3'-P 3 P 4 P 5 C 

when they are hybridized to the target, normalized to the core 
probe. In alternative embodiments, average curves are plotted for 
substitution of all the possible nucleotides at each position (the 
45 "families" of mismatched probes) , or the highest intensity is plotted 
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for each position. Thus, the 0 position on the X axis of the graph 
in Fig. 12B represents no substitution and shows the fluorescence 
intensity due to target hybridization to core probe S-3 ' ~ p 3 p 4 p s p 6 • 
Because all values in Fig. 12B are normalized with respect to this 
5 value, the "no substitution" case has a normalized intensity of l. 

When C is substituted at the 3, 4, 5, and 6 positions, the relative 
intensity values are normally less, because none of these sequences 
are exactly complementary to the target in this example. 

The relative fluorescence intensity of a probe/target 

10 complex with a mismatch at the 3'- or 5' -end is typically higher 

than complexes with mismatches in the center of the probe/target 
complex, because mismatches at the end of the probe tend to be less 
destabilizing than mismatches at the center of the probe/target 
complex. Probe/target complexes with mismatches at the 3 '-end of 

15 the probes may impact hybridization less (and thus have a higher 

fluorescence intensity) than those with mismatches at the 5 '-end of 
the probes, presumably due to the proximity of the 3' -end of the 
probe to the substrate surface in this embodiment. Therefore, a 
curve plotting a normalized factor related to binding affinity versus 

20 mismatch position, tends to have the shape of a "crooked smile," as 

shown in Fig. 12B. 

Using this methodology, one can extend a core sequence by 
examining probes on the array that have the same sequence as the core 
probe except for having been extended at one end and optionally 

25 shortened at the other. These probes are evaluated as candidate 

second core sequences to determine which probes are perfectly 
hybridized to the target. By repetition of this process, one can 
determine the complete nucleotide sequence of the target. 

To illustate the method, Fig. 12C shows the 4 possible, 

30 4-member "left extensions" of the core probe S-3 ' -P 3 P 4 P 5 P 6 . 

As shown, the nucleotide adjacent to the sequence of the target 
complementary to S-3 ' -P 3 P 4 P 5 P 6 is either A, T, C, or G, or there is 
no adjacent nucleotide on the target (i.e., P 3 is complementary to 
the 5' -end of the target). Therefore, the possible left extensions 

35 of the P3P4P5P6 core probe are probes S-3 ' -AP 3 P 4 P 5 , S-3 ' -TP 3 P 4 P 5 , 

S-3' -CP 3 P 4 P 5 , and S-3 ' -GP 3 P 4 P 5 . For the purposes of this 
illustration, T is assumed to be actually "correct," as A is in the 
complementary position in the target nucleic acid. 

The upper left hand plot in Fig. 12D illustrates 

4 0 predicted hybridization data for the mismatch profile of the 

S-3'-AP 3 P 4 P 5 probe, with all data normalized to S-3 ' -AP 3 P 4 P 5 . Data 
points for all substitutions at each of the 2-5 positions are shown, 
but the average data for the three substitutions at each position 
could also be utilized, a single substitution at each position can be 

45 utilized, the highest of the three values may be utilized, or some 
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other combination. As shown in the S-3'-AP 3 P 4 P 5 graph, one point 
shows much higher binding affinity than the rest. This is the T 
substitution for A at position two. The remaining data in the 
AP 3 P 4 P 5 graph have the normal "smile" characteristics shown in 
5 Fig. 12B. Similar plots are developed for the C and G substitutions 

shown in the bottom portion of Fig. 12B. In each case, all 
datapoints are normalized to the presumed "core" probe in the graph. 

The T extension graph, shown in the upper right hand 
portion of Fig. 12D, will not have aberrant curves like the 

10 3'-AP 3 P 4 P 5 graph and others, because none of the monosubstitutions at 

position 2 of the 3'-TP 3 P 4 P 5 probe will be exactly complementary to 
the target. Accordingly, substitutions of A, C, and G at position 2 
all produce the characteristic "smile" plots predicted for probes 
with single base mismatches relative to the target. In addition, the 

15 fluorescence intensity of the T substituted probe/target complex will 

normally be higher than the fluorescence intensity of the C, G, and A 
probe/target complexes . These data can be used in various 
combinations to determine which of the extensions is "correct" 
and thereby determine the sequence of the target nucleic acid. 

20 From the data shown in Figs. 12A to 12D, one concludes 

that the probe exactly complementary to the left extension of the 
target relative to the core probe complementary sequence has an A 
monomer at position 2 in the target. 

This process is repeated until none of the graphs have 

25 appropriate characteristics, at which time it is concluded that an 

end of the target has been reached. Similarly, right extensions are 
evaluated until the end of the target (or end of the sequence of 
interest) is reached. 

The above techniques can obviously be conducted through 

30 manual observation of the hybridization data. However, in preferred 

embodiments the data are analyzed using one or more appropriately 
programmed digital computers. An exemplary system is illustrated 
in Fig. 13. As shown therein, the system includes a computer or 
computers 302 operated under the control of a CPU and including 

35 memory 304, such as a hard disk, and memory 306, such as dynamic 

random access memory. The computer is used to control a scanning 
device 308 that measures the fluorescence intensity or other related 
information from a labeled target nucleotide coupled to portions of a 
substrate 2. The substrate 2 contains probe nucleotides of known 

4 0 sequence at known locations thereon. A user provides input via input 

devices 313. 

Fluorescence intensity or other related information is 
stored in the memory 304/306. CPU 310 processes the fluorescence 
data to provide output to one or both of print device 312 or display 
45 314. The data are processed according to the methods described 
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herein, and output in the form of graphs such as those shown above, 
or in sequence of nucleic acid monomers, or in simple (+)/(-) output, 
or other results of the analysis of such data .may be obtained. 
Suitable computers include, for example, an IBM PC or compatible, a 
5 SPARC workstation, or similar device. 

Fig. 14 is a flowchart for a typical computer program 
used to evaluate an array of n-mers and identify the sequence of an 
exactly complementary (for mismatch analysis) or a larger k-mer (for 
sequencing or other purposes) . As shown therein, the system first 

10 identifies a core probe at step 402 by, for example, selecting a 

probe having the highest binding affinity of some specified set of 
probes. The present method will often be operational in iterative 
processes, where the highest affinity probe in the array is not 
selected after the first iteration, and in other cases, it may be 

15 worthwhile, for example, to select the one, two, three, or more 

strongest binding probes and perform left and right extensions on 
each, then store and compare this information with other data before 
providing the final output. The results can aid in confirmation of 
the correct sequence. 

20 At step 404 the system identifies all left extensions of 

the core n-mer. At step 406 the system selects the appropriate left 
extension by one or both of: 

- determining which of the left extensions exhibits 
the behavior most consistent with a preset monomer 

25 substitution pattern, and/or 

- selecting the left extension exhibiting the 
highest binding affinity. 

The above selection criteria and others may in some embodiments be 
used in an AND fashion, i.e., both of the criteria must be met or the 

30 system assumes that one has either reached the terminal monomer or 

the system is not performing acceptably. In alternative embodiments, 
one of the criteria may be selected as a primary selection mechanism, 
and the other may be used to provide the user with warnings, 
potentially incorrect selections, or alternate selections. 

35 Thereafter, the system determines if the selection 

criteria have met some minimum standard at step 408. If not, then 
the system assumes that the end of the sequence has been reached at 
step 410. If the selection criteria have been met, then the process 
is repeated beginning at step 404 with the new "core" selected as the 

40 correct extension from the previous core. 

Thereafter, the process is effectively repeated for 
right extensions. At step 412 right extensions are identified. At 
step 414 a preset mismatch profile probe is identified and/or high 
affinity right extension. At step 416 the system determines if the 

45 terminus of the molecule has been reached. If not, then the process 
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is repeated to step 412. If so, then the system assumes that the 
molecule has been sequenced, and the process is terminated with 
appropriate output to a printer or other output device . 

5 Another embodiment of a method for analyzing data using 

mismatch analysis is exemplified in Figure 15a, which presents a 
flowchart for a typical computer program designed to determine the 
sequence of an unknown target fragment. Such methods would normally 
be performed in an appropriately programmed digital computer such as 

10 an IBM PC or equivalent, a Sun workstation, or other similar computer 

system. As shown in step 502, hybridization data (such as the 
fluorescence intensity data as described above, although other data 
obtained through other techniques which are well known to those 
skilled in the art may also be used) are entered into the system. In 

15 some embodiments, the data may be input directly from the 

experimental system, while in others, the data may be collected in a 
separate system. At step 504, the system selects probe sequences 
corresponding to the "best" data in the data set. In embodiments 
which utilize fluorescence intensity data, the best data will 

2 0 typically correspond to those probes which exhibit the highest 

intensity levels. A default percentage of the data set may be 
employed, or optionally, the operator may direct the system to choose 
a specific percentage of the most intense data, with the percentage 
chosen based on such factors as the size of the data set, the 

25 system's computing capacity, and the target fragment size. At step 

506, the system identifies the set of all fragments of a specified 
length which can be derived through combination of the probe 
sequences selected in step 504. The fragments thus identified are 
referred to as "candidate" fragments. 

30 At step 508, the system constructs mismatch analysis 

affinity plots as previously described for perfect match probes for 
each candidate fragment in the set, utilizing the measured 
hybridization data from the nucleic acid probes. At step 510, each 
affinity plot is scored to reflect how well it corresponds to the 

35 expected shape of an affinity plot after accounting for cross 

hybridization of the probe with other parts of the target DNA. 

Once the scoring step 510 has been completed, the program 
totals the scores from all the affinity plots for each candidate 
fragment (step 512) , and proceeds to examine the next candidate 

40 fragment at step 514. At step 516, when all the candidate fragments 

have been examined, the system compares the totalled affinity plot 
scores for the individual candidate fragments. Finally, at step 518, 
the system selects the candidate fragment or fragments for which the 
score indicates the closest overall match with the template. The 



CGK00035958 



WO 95/00530 



PCT/US94/07106 



24 

sequence of this candidate fragment corresponds to the sequence of 
the unknown target fragment. 

Figure 15b is a flowchart for one preferred embodiment of 
a method to score the affinity plots. Initially, a predetermined 
5 plot of expected binding affinity versus mismatch position 

("template") is entered into the system (step 602). Typically, this 
template will exhibit the "smile" shape as previously described. At 
step 604, the system compares each affinity plot with the template as 
specified in steps 606 through 620. In steps 606 and 608 the system 

10 compares the affinity value for each mismatch position with the value 

of the corresponding mismatch position in the template. At step 610, 
the affinity value is examined, initially, to determine whether the 
measured value falls within a chosen interval of the template value. 
In some embodiments, this interval is set to a default level; in 

15 others, it is selected by the operator, based on factors such as the 

expected magnitude of the experimental error. 

If the value at a given mismatch position does fall 
within the interval, the process advances to step 618, in which the 
system marks that mismatch position for inclusion in subsequent 

20 calculations. If the measured value at a mismatch position does not 

fall within the chosen interval, then at step 612, the system 
compares the base sequence including the mismatch ("mismatched 
sequence") with the entire candidate fragment. At step 614, the 
process determines whether there is a perfect match between that 

25 mismatched sequence and any portion of the candidate fragment. If 

such a match is found, the system excludes that mismatch position 
from its scoring of the plot (step 616), and continues to step 620; 
otherwise, that position is marked for inclusion in subsequent 
calculations (step 618) . At step 620, the process continues until 

30 each mismatch position has been examined. The fit of the 

experimental values with the template values is then determined in 
the calculation of step 622 for each mismatch position which had been 
marked for inclusion in step 618. In some embodiments, this 
calculation may include scoring each point as either in the interval 

35 or outside the interval; other embodiments may utilize methods to 

whereby each included position is scored to reflect its divergence 
from the anticipated value, including such well known techniques as 
the Root-Mean-Squares analysis. Finally, the system calculates a 
total score for each affinity plot in step 624. 

40 

Example 

An array of 8-mer probes was constructed according to the 
methods described in U.S. Patent No. 5,143,854 issued to Pirrung et 
45 aJK and in U.S. Application Serial No. 07/805,727, both incorporated 
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herein by reference for all purposes. The array was incubated with a 
16-mer target, and scanned for fluorescence intensity utilizing known 
methods. ( See, e.g. , U.S. Application Serial No. 08/195,889, 
incorporated herein by reference for all purposes) . The resulting 
5 fluorescence intensity data were analyzed according to the foregoing 

method. After the affinity plots had been created and scored, four 
candidate fragments with scores ranging from 69 - 71 were identified: 
5 ' - AGTTGTAGTGGATGGT , TGTTGTAGTGGATGGT , GGTTGTAGTGGATGGT , and 
CGTTGTAGTGGATGGT . The four likely fragments differ only in the 

10 identity of the base at the 5' end of the fragment, which does not 

receive the benefit of the duplicative analyses that the other 
positions receive. The separation of the scores for these four 
candidates from the scores for the other candidate fragments was 
substantial; the next highest scoring candidates received a score of 

15 54 . 

In yet another embodiment of the present invention, the 
arrays discussed above can be used to identify the sequence of an 
unknown target fragment by the following method. 

A set of hybridization data from probes of length "n" 

20 (where "n M is equal to the number of bases in the probe and is less 

than the length of the unknown target fragment) is entered into an 
appropriately programmed computer system. The hybridization data 
sets contemplated herein include sets of fluorescence intensity data 
as previously described, as well sets of data obtained through other 

25 techniques known to those skilled in the art. The system constructs 

a directed graph (as described in Harary, F., Graph Theory . Addison- 
Wesley, Reading, MA (1969) and in Ahuja, R.K., Magnanti, T.L., Orlin, 
J.B., Network Flows, Theory, Algorithms, and Applications . Prentice 
Hall, New York (1993), both incorporated herein for all purposes), 

30 wherein the vertices (or "nodes") are the base sequences 

corresponding to all (n-l)-mers of the data set, and the edges are 
the base sequences corresponding to all n-mers of the set. Each edge 
is constructed such that the edge connects from the initial contained 
(n-l)-mer to the terminal contained (n-l)-mer. A sample of such a 

35 directed graph where the target k-mer is known to be ACTGTTG and the 

n-mers are 3-mers is included as Figure 16. In Figure 16, only the 
edges corresponding to the 3-mer sequences are included; in the 
typical case in which the target base sequence is unknown, the system 
would include all possible edges connecting the n-mer nodes in its 

40 analysis. 

In addition, the system creates a source node in the 
directed graph, distinct from any of the (n-l)-mer nodes, and an edge 
from that source node to each of the (n-l)-mer nodes in the graph. 
Finally, the system creates a sink node in the graph, distinct from 
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any of the previously created nodes, and an edge from each of the (n- 
l)-mer nodes of the graph to that sink node. 

The system then assigns a value to, each of the graph's 
edges {"cost") . The edges leading from the start node and the edges 
5 leading to the sink node are all assigned a cost of zero. All other 

edges are assigned a cost as a function of the probability that the 
base sequence corresponding to that edge is a perfect match of a 
segment in the target k-mer. This probability is referred to as the 
"probability of a perfect match". In a preferred embodiment, the 

10 cost assigned to each n-mer edge is equal to the negative of the 

natural log of the probability that the sequence corresponding to 
that edge is a perfect match of a segment in the target fragment 
divided by the probability that the sequence corresponding to that 
edge is not a perfect match. In some embodiments, the edges may be 

15 assigned a separate, additional value ("capacity"); for such 

embodiments, the edges leading from the start node to each of the 
other nodes and the edges leading to the sink node are all assigned a 
capacity of one; all other edges of the graph are assigned infinite 
capacity. 

20 Beginning at the start node and following along the 

directional edges, numerous paths can be traced through the graph, 
ending at the sink node, which may possibly correspond to the 
sequence of the target fragment . Since each edge of the graph has 
been assigned a cost, the total cost for each of these paths can be 

25 determined. In the final step of this method, the system applies an 

algorithm to the graph to determine this total cost for each possible 
path through the graph, beginning at the start node and ending at the 
sink node. Algorithms contemplated for use in this step include 
minimum cost /maximum flow programming algorithms well known to those 

30 skilled in the art, such as described in Ahuja et al . , Network Flows. 

Theory, Algorithms, and Applications , Prentice Hall, New York (1993), 
incorporated herein by reference for all purposes. Using such an 
algorithm, the system selects the path with total cost for which 
there is the highest combined probability of perfect matches. For 

35 example, in embodiments utilizing a minimum cost /maximum flow 

algorithm, the path with the highest probability of perfect match 
corresponds to the path with the minimum cost. The sequence of bases 
associated with that path, therefore, corresponds to the most likely 
sequence for the target . 

40 In a preferred embodiment, the probability of a perfect 

match for an individual probe is determined from fluorescence 
hybridization data as follows. For each probe of the chosen array, 
the hybridization data for a set of related probes are entered into 
the system; the pattern of intensities of the probes related to the 

45 probe of interest is then analyzed for match with the known target 
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sequence utilizing the mismatch analysis described above. The set of 
related probes contemplated by this step would typically include 
those probes with one or two base mismatches, as chosen by the 
operator based upon factors such as computing capacity and the size 
5 of the data set available. The data included in the set would 

typically include the intensities and the base composition of the 
probes, and may optionally include data in regard to base 
substitutions which would cause the probe to become a perfect match 
elsewhere in the target ("crosstalk") . Preferably, these data will 

10 be collected for a large series of experiments utilizing known DNA 

sequences as the target; the merged data from the different 
experiments can then be analyzed through a Neural Network, or other 
such well known learning based classification method (as described 
in, for example, Rich, E . , Knight, K. Artificial Intelligence . 2nd 

15 ed., McGraw-Hill, Inc., New York (1991), which is incorporated herein 

for all purposes) , to develop an algorithm to predict whether a probe 
is a perfect match as a function of intensity, base composition, and, 
optionally, crosstalk. The confidence level of such a prediction can 
be estimated in the algorithm; this confidence level corresponds to 

20 the probability of a perfect match that is used in assigning cost 

values to the edges of the directed graph. 



D . Applications 

25 The techniques described herein will have a wide range of 

applications, particularly wherever desired to determine if a target 
nucleic acid has a particular nucleotide sequence or some other 
sequence differing from a known sequence. For example, one 
application of the inventions herein is found in mutation detection. 

30 These techniques may be applied in a wide variety of fields including 

diagnostics, forensics, bioanalytics , and others. 

For example, assume a "wild- type" nucleic acid has the 
sequence 5'-N 1 N 2 N 3 N 4 where, again, N refers to a monomer such as a 
nucleotide in a nucleic acid and the subscript refers to position 

35 number. Assume that a target nucleic acid is to be evaluated to 

determine if it is the same as 5'-N 1 N 2 N 3 N 4 or if it differs from this 
sequence, and so contains a mutation or mutant sequence. The target 
nucleic acid is initially exposed to an array of typically shorter 
probes, as discussed above. Thereafter, one or more "core" sequences 

4 0 are identified, each of which would be expected to have a high 

binding affinity to the target, if the target does not contain a 
mutant sequence or mutation. In this particular example, one probe 
that would be expected to exhibit high binding affinity would be the 
complement to 5'-N 1 N 2 N3 (3'-P 1 P 2 P 3 ), assuming a 3-mer array is 
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utilized. Again, it will be recognized that the probes and/or the 
target may be part of a longer nucleic acid molecule. 

As an initial screening tool, the absolute binding 
affinity of the target to the 3'-P 1 P 2 P 3 probe will be utilized to 
determine if the first three positions of the target are of the 
expected sequence. If the complement to 5'-N 1 N 2 N 3 does not exhibit 
strong binding to the target, it can be properly concluded that the 
target is not of the wild-type. 

The single base mismatch profile can also be utilized 
according to the present invention to determine if the target 
contains a mutant or wild- type sequence. Figs. 17A and 17B 
illustrate typical illustrative plots resulting from targets that are 
wild- type (Fig. 17A) and mutant {Fig. 17B) . As shown, the single 
base mismatch plots for wild-type targets generally follow the 
15 typical, smile-shaped plot. Conversely, when the target has a 

mutation at a particular position, not only will the absolute binding 
affinity of the target to a particular core probe be less, but the 
single base mismatch characteristics will deviate from expected 
behavior . 

20 According to one aspect of the invention, a substrate 

having a selected group of nucleic acids (otherwise referred to 
herein as a "library" of nucleic acids") is used in the determination 
of whether a particular nucleic acid is the same or different than a 
wild- type or other expected nucleic acid. Libraries of nucleic acids 

25 will normally be provided as an array of probes or "probe array." 

Such probe arrays are preferably formed on a single substrate in 
which the identity of a probe is determined by ways of its location 
on the substrate. Optionally, such substrates will not only 
determine if the nucleotide sequence of a target is the same as the 

30 wild-type, but it will also provide sequence information regarding 

the target. Such substrates will find use in fields noted above such 
as in forensics, diagnostics, and others. Merely by way of specific 
example, the invention may be utilized in diagnostics associated with 
sickle cell anemia detection, detection of any of the large number 

35 of P-53 mutations, for any of the large number of cystic fibrosis 

mutations, for any particular variant sequence associated with the 
highly polymorphic HLA class 1 or class 2 genes (particularly class 2 
DP, DQ and DR beta genes), as well as many other sequences associated 
with genetic diseases, genetic predisposition, and genetic 

4 0 evaluation. 

When a substrate is to be used in such applications, it 
is not necessary to provide all of the possible nucleic acids of a 
particular length on the substrate. Instead, it will be necessary 
using the present invention to provide only a relatively small subset 
45 of all the possible sequences. For example, suppose a target nucleic 
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acid comprises a 5-base sequence of particular interest and that one 
wishes to develop a substrate that may be used to detect a single 
substitution in the 5-base sequence. According to one aspect of the 
invention, the substrate will be formed with the expected 5-base 
5 sequence formed on a surface thereof, along with all or most of the 

single base mismatch probes of the 5-base sequence. Accordingly, 
it will not be necessary to include all possible 5-base sequences 
on the substrate, although larger arrays will often be preferred. 
Typically, the length of the nucleic acid probes on the substrate 

10 according to the present invention will be between about 5 and 100 

bases, between about 5 and 50 bases, between about 8 and 30 bases, or 
between about 8 and 15 bases. 

By selection of the single base mismatch probes among all 
possible probes of a certain length, the number of probes on the 

15 substrate can be greatly limited. For example, in a 3 -base sequence 

there are 69 possible DNA base sequences, but there will be only one 
exact complement to an expected sequence and 9 possible single base 
mismatch probes. By selecting only these probes, the diversity 
necessary for screening will be reduced. Preferably, but not 

20 necessarily, all of such single base mismatch probes are synthesized 

on a single substrate. While substrates will often be formed 
including other probes of interest in addition to the single base 
mismatches, such substrates will normally still have less than 50% of 
all the possible probes of n-bases, often less than 30% of all the 

25 possible probes of n-bases, often less than 20% of all the possible 

probes of n-bases, often less than 10% of the possible probes of 
n-bases, and often less than 5% of the possible probes of n-bases. 

Nucleic acid probes will often be provided in a kit for 
analysis of a specific genetic sequence. According to one embodiment 

30 the kits will include a probe complementary to a target nucleic acid 

of interest. In addition, the kit will include single base 
mismatches of the target. The kit will normally include one or more 
of C, G, T, A and/or U single base mismatches of such probe. Such 
kits will often be provided with appropriate instructions for use of 

35 the complementary probe and single base mismatches in determining the 

sequence of a particular nucleic acid sample in accordance with the 
teachings herein. According to one aspect of the invention, the kit 
provides for the complement to the target, along with only the single 
base mismatches. Such kits will often be utilized in assessing a 

4 0 particular sample of genetic material to determine if it indicates a 

particular genetic characteristic. For example, such kits may be 
utilized in the evaluation of a sample as mentioned above in the 
detection of sickle cell anemia, detection of any of the large number 
of P-53 mutations, detection of the large number of cystic fibrosis 

45 mutations, detection of particular variant sequence associated with 
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the highly polymorphic HLA class 1 or class 2 genes (particularly 
class 2 DP, DQ and DR beta genes) , as well as detection of many other 
sequences associated with genetic diseases, genetic predisposition, 
and genetic evaluation. 
5 Accordingly, it is seen that substrates with probes 

selected according to the present invention will be capable of 
performing many mutation detection and other functions, but will need 
only a limited number of probes to perform such functions. 

10 Examples 

1 . (G+T)** Array and Differential Sequencing 

A (G+T) 8 array was prepared and incubated with 
1 nM 5' -AACCCAACCCC- fluorescein (representing a mutant sequence when 
compared to 5 ' -AACCCAAACCC) , and scanned to test whether the sequence 

15 was "wild" or "mutant." The resulting image is given in Fig. 18. 

Four overlapping, exactly complementary octanucleotide probe/target 
hybridizations are expected if one is assuming the target should be 
5 ' -AACCCAAACCC with probes: S-3 ' -TTGGGTTG, TGGGTTGG, GGGTTGGG , and 
GGTTGGGG . The results demonstrated that the effect of a single base 

20 change is quite dramatic, especially in the number and identity of 

the different mismatched probe/target complexes that form on the 
array. If one assumes the target nucleic acid generating the signal 
in Fig. 18 is 5 ' -AACCCAAACCC , (i.e., the wild- type) then the mismatch 
profiles for the complementary probe S - 3 ' - TTGGGTTT are shown in 

25 Fig. 19A. The mismatch profile does not have the expected shape, and 

the probe/target complex has a low fluorescence intensity. The 
strong peak corresponding to a mismatch in position 8 indicates that 
the "correct" base in this position in the target is probably an A, 
because only A and C are found in the target in this experiment. 

30 Mismatch position 6 also shows a small peak. By contrast, a similar 

plot using the probe sequence S - 3 TTGGGTTG probe sequence as a core 
yielded the "smile" shape and high fluorescence intensity. In 
Fig. 19B the same profile for the next 8-mer probe is shown. The 
peaks have shifted one position to the left, again confirming that 

35 the sequence varies from wild-type at position 8 in the target. 

These correspond to the same positions in the original 11-mer target 
fragment. These data predict that there is a single base change in 
position 8 of the target, as compared to the wild- type. 

All of the mismatch probe profiles corresponding to the 

40 assumed fragment 5 ' -AACCCAACCCC, are shown in Fig. 19C. One observes 

the mutant position "moving" down the sequence. Finally, in Fig. 19D 
the mismatch plots are shown corresponding to the four probes that 
complement 5 ' -AACCCAACCCC , with the expected smile characteristics. 
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E . Conclusion 

The present inventions provide improved methods and 
devices for the study of nucleotide sequences and nucleic acid 
interactions with other molecules. The above description is 
illustrative and not restrictive. Many variations of the invention 
will become apparent to those of skill in the art upon review of thi 
disclosure. Merely by way of example certain of the inventions 
described herein will have application to other polymers such as 
peptides and proteins, and can utilize other synthesis techniques. 
The scope of the invention should, therefore, be determined not with 
reference to the above description, but instead should be determined 
with reference to the appended claims along with their full scope of 
equivalents . 



20 
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WHAT IS CLAIMED IS: 

1. A method of sequencing a target nucleic acid with a 
plurality of nucleic acid probes, said probes having fewer bases than 

5 said target, comprising the steps of: 

contacting said probes with said targets- 
identifying a first probe that specifically hybridizes to 

said targets- 
selecting a first set of extension probes that comprise 
10 at least two of A, C, T, U, and G extensions of said first probe; and 

identifying one of said first set of extension probes 
that hybridizes specifically to said target more strongly than others 
of said first set of extension probes, whereby said one of said 
extension probes identifies a base in said target nucleic acid. 

15 

2. A method as recited in claim l wherein substantially all 
of said nucleic acid probes comprise n nucleotides, and wherein said 
extension probes comprise n-1 nucleotides of said first probe . 

20 3. A method as recited in claim 1 further comprising the 

steps Of : 

selecting a second set of A, C, T, U, and G extension 
probes that extend in a direction opposite of said first set of 
extension probes; and 
25 identifying one of said second set of extension probes 

that hybridizes specifically to said target more strongly than others 
of said second set of extension probes, whereby said one of said 
second set of extension probes identifies a second base in said 
target nucleic acid. 



30 



35 



4. The method as recited in claim 1 further comprising the 
step of : 

repeating said steps of selecting sets of extension 
probes and identifying extension probes five or more times. 



5. The method as recited in claim 1 wherein said step of 
identifying further comprises the steps of: 

identifying single base mismatch probes, said single 
base mismatch probes comprising at least two of A, C, T, U, and G 
4 0 monosubstitutions of said first set of extension probes; 

recording hybridization affinity data of said single base 
mismatch probes; and 

selecting one of said first set of extension probes 
as a correct extension of said first probe when said hybridization 
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affinity data conform to expected hybridization affinity data of said 
single base mismatch probes. 

6. The method as recited in claim 5 wherein said expected 
hybridization data comprise: 

higher binding affinity for probe/target complexes with a 
mismatch at termini of said extension probes; and 

lower binding affinity for probe/target complexes with a 
mismatch at internal portions of said complexes. 



7. The method as recited in claim 6 wherein: 

said hybridization data are normalized to a hybridization 
value for one of said extension probes; and 

said step of identifying comprises selecting one of 
15 said extension probes having terminal single base mismatch probes 

that do not have normalized hybridization values higher than a 
normalized value of said one of said extension probes. 

8 . The method as recited in claim 1 wherein the step of 
20 identifying comprises the step of selecting one of said set of 

extension probes that exhibits a higher binding affinity to said 
target than other extension probes . 

9 . The method as recited in claim 1 wherein said step of 
25 identifying is conducted in an appropriately programmed computer. 

10. A method of determining if a nucleotide sequence of a 
target nucleic acid is the same as a sequence of a first nucleic acid 
comprising : 

30 contacting said target nucleic acid to a plurality of 

nucleic acid probes; 

determining the affinity of said target to probes 
identical to, but for a single base mismatch, of said subsequence; 
and 

35 determining that said nucleotide sequence of said target 

is the same as said first nucleic acid if said affinity of said 
target to probes identical to but for a single base mismatch follows 
a predetermined pattern. 

4 0 11. The method as recited in claim 10 wherein said 

predetermined pattern comprises affinity of said single base mismatch 
probes normalized to affinity of a perfect complement of said 
subsequence. 
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12. The method as recited in claim 11 wherein said affinity 
of single base mismatch probes are plotted as affinity versus 
mismatch position, and normalized to said affinity of a perfect 
complement of said subsequence. 

5 

13. The method as recited in claim 10 further comprising the 
step of determing that said nucleotide sequence of said target is not 
the same as said first nucleic acid if said affinity of said target 
to probes complementary to single base mismatches does not follow a 

10 predetermined pattern. 

14. A probe array of nucleic acids, said probe array selected 
from all possible probes to comprise an exact complement to a target 
nucleic acid, and single base mismatches of said exact complement. 

15 

15. A library as recited in claim 14 wherein said nucleic 
acid probes are of a length between about 8 and 15 bases. 

16. A library as recited in claim 14 wherein said library is 
20 on a single substrate. 

17. A library as recited in claim 14 wherein said library 
comprises probes of n-bases or less, and wherein said library 
comprises less than 50% of all possible probes of n-bases. 

25 

18. A library as recited in claim 14 wherein said library 
comprises probes of n-bases or less, and wherein said library 
comprises less than 10% of all possible probes of n-bases. 

30 19. A nucleic acid probe kit comprising a core nucleic acid 

probe, said core probe exactly complementary to a nucleic acid 
target, and selected A, C, T, U, and G single base substitutions of 
said core probe . 

35 20. A nucleic acid probe kit as recited in claim 19 consisting 

essentially of said core probe and A, C, T, and G single base 
substitutions of said core probe. 

21. A nucleic acid probe kit as recited in claim 19 further 
40 comprising instructions for determining if a target sample is the 

same as or different than said target. 

22. A nucleic acid probe kit as recited in claim 19 wherein 
said core probe comprises between 8 and 15 bases. 

45 
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23. A nucleic acid probe kit as recited in claim 19 wherein 
said probes are selected to evaluate a target sample for a genetic 
characteristic selected from the group consisting of sickle cell 
anemia, P-53 mutations, cystic fibrosis mutations, HLA class 1 genes, 
5 and HLA class 2 genes . 



24 . A nucleic acid probe kit as recited in claim 19 wherein 
said probes are selected to evaluate a target sample for sickle cell 
anemia . 



10 



25. A method of sequencing a target nucleic acid comprising 
the steps of : 

contacting said target with an array of probes ; 
identifying selected high affinity probes in said array 
15 with at least an identified hybridization affinity level; 

identifying a plurality of candidate target sequences 
derivable from said high affinity probes; 

identifying mismatch probes wherein said mismatch probes 
contain identical base sequence as said high affinity probe but for 
20 at least one base mismatches; 

comparing an affinity pattern for said mismatch probes to 
an expected mismatch affinity pattern for said candidate target 
sequences ; and 

selecting a target sequence from said candidate target 
25 sequences, said selected target sequence having a best mismatch 

pattern . 

26. The method as recited in claim 25 wherein said expected 
mismatch affinity pattern comprises: 

30 higher binding affinity for probe/target complexes with a 

mismatch at termini of said extension probes; and 

lower binding affinity for probe/target complexes with a 
mismatch at internal portions of said complexes. 

35 27. The method as recited in claim 25 wherein said step of 

comparing affinity pattern comprises the steps of: 

identifying an expected affinity value for each mismatch 
position in said expected mismatch affinity pattern; and 

selecting a plurality of said mismatch positions having 
4 0 an affinity value no more than an identified variation from said 

expected affinity value. 

28. A method of sequencing a target nucleic acid comprising 
the steps of : 

45 contacting said target with an array of shorter probes; 
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36 

determining probability of a perfect match with a portion 
of said target for each of said probes; 

determining a total of said probabilities of a perfect 
match for all combinations of said probes; and 

selecting an identified combination of said probes from 
said combinations of said probes for which said total of 
probabilities of a perfect match is maximized, whereby said 
identified combination of said probes identifies sequence of said 
target nucleic acid. 



29. The method as recited in claim 28 wherein said step of 
determining probability of a perfect match comprises the steps of: 
collecting hybridization data for a plurality of known 

fragments ; 

15 identifying algorithm to predict known probability of a 

perfect match from said hybridization data; and 

identifying said probability of a perfect match for said 
probes using said algorithm. 

20 



25 



30 



35 
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